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Abstract 

We compute several matrix elements of dimension-six four-fermion operators and 
extract their l?-parameters. The calculations have been performed with the tree-level 
Clover action at (3 = 6.0. The renormalization constants and mixing coefficients of the 
lattice operators have been obtained non-perturbatively. In the MS renormalization 
scheme, at a renormalization scale fj, ~ 2 GeV, we find Bk {B^ 2 ) = 0.66(11), B^ 2 = 

0.72(5) and B^ 2 = 1.03(3). The result for B^ 2 has important implications for the 
calculation of e'/e. 



1 Introduction 



The lattice evaluation of matrix elements of operators between hadronic states is a 
necessary ingredient to the calculation, from first principles, of a wide class of physical 
observables. In the effective Hamiltonian approach, weak amplitudes are expressed in 
terms of perturbative Wilson coefficients multiplied by matrix elements of local opera- 
tors, which can be evaluated on the lattice. The results are conventionally presented in 
terms of i?-parameters, which measure the deviation of the matrix elements from their 
value in the Vacuum Saturation Approximation (VSA). These quantities are subject 
to significant QCD corrections. 

In this paper we focus on AS = 2 and AI = 3/2 transition amplitudes. The 
former case is characteristic of K°—K° oscillations, which are related to indirect CP 
violation, parametrized by e. We compute Bk, the ^-parameter of the matrix element 
(K \O^ s=2 \K°) , where the dimension-six, four-fermion operator A5=2 has a "left- 
left" chiral structure (see sec. |2|). We also compute the AI = 3/2 amplitudes relevant 
in K — > 7T7T decays. These amplitudes are essential for theoretical predictions of direct 
CP-violation, parametrized by e'. Matrix elements of the form (7nr\0 3 / 2 \K) enter 
in the calculation of e', with two of the electro-penguin operators (O^ 2 and O^ 2 ) 

having a "left-right" chiral structure, and one operator (0 9 ) having a "left-left" one. 
Using Chiral Perturbation Theory (CPTh), these matrix elements can be related to 
the single-state matrix elements (tt + \0 3 ^ 2 \K + ) . In this work we compute the latter, 
parametrized in terms of l?-parameters denoted as B^ 2 and B^ 2 (the conventional 
basis of operators for the AS* = 1 effective Hamiltonian can be found in refs. jj]; we 
also remind the reader that B 9 = Bk). 

The Wilson lattice regularization breaks chiral symmetry. This implies that the 
AS = 2 operator, which would otherwise renormalize multiplicatively, mixes with 
operators belonging to different chiral representations |^, ||. The same is true for 
the two AI = 3/2 operators, which would otherwise only mix with each other, as 
they belong to the same chiral representation. Because of the mixing induced by the 
lattice, the correct chiral behaviour of the operators is achieved with Wilson fermions 
only in the continuum limit. For the AS = 2 case, for example, restoration of the 
chiral properties amounts to the vanishing of the matrix element as txik — ► 0. In 
practice, the mixing of the AS = 2 operator with operators of "wrong" naive chirality, 
computed at small but finite cutoff a -1 , spoils the expected chiral behaviour [|], [|] []. 
This problem can be attributed to two sources of systematic error in the computation 
of the matrix elements, namely the determination of the mixing coefficients in one- 
loop perturbation theory, and the 0(a) (a is the lattice spacing) discretization errors. 
Several methods have been proposed in order to improve the determination of the 
mixing coefficients. One of them consists in evaluating the renormalization constants 

1 In the Staggered fermion approach, where chiral symmetry is partially preserved, the AS = 2 matrix 
element displays the correct chiral behaviour. Thus, the -Bjf-parameter obtained with staggered fermions Q 
has been deemed more reliable. 
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computed in Standard Perturbation Theory (SPT) using an effective coupling Q which 
should reduce higher order corrections; we refer to it as Boosted Perturbation Theory 
(BPT). Another is the Non-Perturbative Method (NPM) for the computation of the 
renormalization constants on quark and gluon external states, as proposed in ref. ||. 
Finally, in the spirit of ref. §, the lattice mixing coefficients can also be obtained 
non-perturbatively by using the Ward Identity Method (WI's), with external quark 
and gluon states ||; see also refs. [ 10 1 . On the other hand, recent studies [Il]-[|l3| (see 
also the present work) have shown that reducing the discretization error by using the 
tree-level Clover action does not improve the chiral behaviour of the matrix elements, 
even when BPT is implemented in the definition of the renormalized AS" = 2 operator. 
Instead, a good chiral behaviour has been observed by evaluating the renormalization 
constants non-perturbatively, either with the NPM (refs. [|ll[-[13] and this work) or 
by using WI's (ref. |||). The fact that the restoration of the correct chiral behaviour 
has been seen both with the tree-level Clover action (refs. Ji I Jl^] and the present 
work) and with the Wilson action (ref. ||), suggests that discretization effects are less 
important than those due to the perturbative evaluation of the mixing coefficients. 
Therefore, the recent Wilson fermion estimates of Bk from the NPM or the WI's, 
are considerably more reliable than those of earlier studies, based on perturbative 
renormalization. 

Given the success of the NPM in the computation of Bk, we also apply it to the 
evaluation of the two AI = 3/2 ^-parameters of the electro-penguin operators, B 7 

Q /9 

and B s . A recent calculation, using the Wilson action and with the renormalization 
constants obtained in BPT, found B^ 2 = 0.81(3) fll4j] , as opposed to the earlier re- 
sults B^ 2 ~ 1 of refs. [15, 16]. In the present work, using the Clover action and the 
NPM for the evaluation of the renormalization constants we find B£ = 1.03(3). No- 
tice that, although obtained with a different lattice fermion action, our BPT estimate 
-Bg = 0.83(2) is, instead, fully compatible with that of ref. [1J]. Our preferred value 
B s = 1.03(3), obtained with an improved operator renormalized non-perturbatively, 
only suffers from discretization errors which are O(gQa). We thus believe that it is 
more reliable than previous results which have been obtained with the (non-improved) 
operator renormalized in one-loop perturbation theory. The same situation character- 

3 II 

izes B 7 : our NPM and BPT results are not in agreement, but the latter is compatible 
with the value quoted in |14[] . We stress that a precise determination of B s , com- 
bined with an equally reliable estimate of the strange quark mass, is essential to the 
determination of the ratio e'/e. For example, the uncertainties in the measurement 
of B^ 2 , combined with the controversial results for m s |17]-[18| may change e'/e up 

3 /2 

to a factor of 2 to 3. In view of its importance, we believe that B 8 should also be 
computed by applying the NPM with the Wilson action also and/or by using WI's for 
the determination of the lattice mixing coefficients of the renormalized operator. 

An extensive study of the renormalization properties of the four-fermion operators 
can be found in [IIJ] . There we detail all the theoretical and numerical issues of relevance 
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to the non-perturbative renormalization of the AS = 2 and AI = 3/2 operators. We 
have used these results in the present study. 

The paper is organized as follows: in sec. ^ we introduce the operators of inter- 
est and review briefly their phenomenological implications; in sec. ||] we address the 
problem of operator mixing and give a brief account of the NPM for the computation 
of the renormalization constants; in sec. |] define the i?-parameters and discuss their 
extraction from lattice correlation functions; in sec. ||, we present our results for Br, 
B 3 ^ 2 and B^ 2 ; finally, in sec. |6| we compare our results to those previously obtained 
in the literature. 



2 AS = 2 and AI = 3/2 transitions 

Schematically, the AS = 2 effective Hamiltonian has the form 



H AS = 2 = ^M" w 6 AS =^Mx c , x t , M ) (1) 



where x q = m q jM^ and $ is a combination of the Inami-Lim functions [20] weighted 
by the CKM matrix elements. All perturbative QCD corrections (known at the next- 
to-leading order (NLO) ]21|| ) are included in (which is, hence, //-dependent). The 
rest of the notation is standard: Gf is the weak Fermi coupling, Mw the mass of the W 
boson and /i the renormalization scale. The AS = 2 operator in the above expression 
is the renormalized one. It is defined as follows: 

O as = 2 = srfdsrfd, (2) 

where 7^ = 7 M (1 — 75). The CP-violation parameter €k is defined as: 

where Amx is the K®-Kg mass splitting, mx and fx denote the if-meson mass and 
decay constant, respectively. 

In the AS* = 1 case, the AI = 3/2 contribution to e' can be written as 

W 3 / 2 oc -G F lm\ t [c 7 (0 3 7 /2 )vsaB 3 7 /2 + C s (0 3 8 /2 ) V saB 3 8 /2 + C 9 (0 3 9 /2 ) VS aB 3 9 /2 ] , 

(4) 

where At = V t * s Vtd contains the CKM matrix dependence and {0 3 ^ 2 )vsa stands for the 
VSA matrix element of the corresponding operator. The definitions of the operators 
can be found for example in refs. [^, p3j ]. The Wilson coefficients, known to NLO [22] 
[25], are denoted by Ck = Ck(M\y/ /i) and the operators Ok (k = 7,8,9) are defined 
as: 

3/2 = (s a ^d a ){u^Uj3-d^d } + (s a ^u a )(up^dp), 

3 8 /2 = (s a ^dp){up^u a -d^d a } + (s a ^u p )(u ^d a ), (5) 

Ol /2 = (s a rfd a ){uprfuf3 - dprfdp} + (s a rfu a )(uprfdp), 
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where 7^ = 7^(1 + 75) and a, (3 = 1-3 are colour indices. The definitions of the 
.B-parameters of eqs. @ and (||) will be given in sec. |j. 

Two observations are necessary at this point. Firstly, since we are interested in 
computing the matrix elements {K°\6 AS=2 \K°) and (tt + \61 /2 \K + ) (with k = 7,8,9), 
only the parity-conserving parts of the operators of eqs. @ and (||) enter in the calcu- 
lation. Secondly, on the lattice the above matrix elements are obtained in the standard 
way by studying the asymptotic behaviour, at large time separations, of hadronic cor- 
relation functions of the form (P(y)O(0)P(x)) (see eqs. ([14]) below), with P denoting 
suitable pseudoscalar densities which we use as meson sources and sinks. The Wick 
contractions of the quark fields in the correlation functions give rise to diagrams which 
are both "eight" -shaped and "eye"-shaped. The latter, however, cancel in the limit of 
degenerate up and down quark masses. Since our results are obtained in this limit, 
complicated subtractions of lower dimensional operators, necessary for the removal of 
the power divergences of the "eye" -diagrams, are avoided. 



3 Non-Perturbative renormalization 

The NPM for the evaluation of the renormalization constants of lattice operators con- 
sists in imposing suitable renormalization conditions on lattice amputated quark cor- 
relation functions ||. In our case, we compute four-fermion Green functions in the 
Landau gauge. All external quark lines are at equal momentum p. After amputating 



and projecting these correlation functions (see refs. and [19] for details), the renor- 
malization conditions are imposed in the deep Euclidean region at the scale p 2 = fj, 2 . 
This renormalization scheme has been recently called the Regularization Independent 
(RI) scheme |2^] (MOM in the early literature) in order to emphasize that the renormal- 
ization conditions are independent of the regularization scheme, although they depend 
on the external states used in the renormalization procedure (and on the gauge). Thus, 
at fixed cutoff (i.e. fixed (3), we compute non-perturbatively the renormalization con- 
stants and the renormalized operator O r1 '(//) in the RI scheme. In order to obtain 
the physical amplitudes, which are renormalization group invariant and scheme inde- 
pendent, the renormalized matrix elements must subsequently be combined with the 
corresponding Wilson coefficients of the effective Hamiltonian. The latter are known 



in continuum perturbation theory, at the NLO, both in the MS scheme ]21j-[|2j|] and 
in the RI scheme p2[ . Although not strictly necessary, since the standard practice 
consists in giving the £>-parameters in the MS scheme, we will express our results both 
in the MS_and RI schemes. In order to obtain the corresponding operators in the MS 
scheme, O ms (/x), the matrix elements of the operators O ri {ix) must be corrected by 



finite matching coefficients [22|. We stress that, if fi is much larger than ^qcd-, the 
NPM of ref. [11] and the WI method of ref. ||, used for the computation of the mixing 
coefficients of the lattice operators, are equivalent, in the chiral limit. This has been 
shown for two-fermion operators in ref. ||; for four-fermion operators it is discussed in 
detail in refs. ||, || (see also ref. §). 
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In |L9[, we have determined non-perturbatively the operator mixing for the complete 
basis of four-fermion operators, with the aid of the discrete symmetries (parity, charge 
conjugation and switching of flavours), in the spirit of ref. Q. The renormalization 
of the parity-conserving operators, relevant to this work, is conveniently expressed in 
terms of the following basis of five operators: 

Q 1 = V xV + Ax A, 
Q 2 = V xV - Ax A, 

Q 3 = SxS-PxP, (6) 
Q 4 = S xS + Px P, 
Q 5 = TxT. 

The operators Qi, ■ ■ ■ ,Q§ form a complete basis on the lattice. In these expressions, 
r x T (with r = V, A, S, P, T a generic Dirac matrix) stands for ^('ipiFipz'ipsT^ + 
^iTifj^ip^T^), where fy, i = 1, ... ,4 are fermion fields with flavours chosen so as to 



(7) 



reproduce the desired operators (see ref. [19] for details). More specifically, the parity- 
conserving component of the four-fermion operator 0^ s=2 corresponds to Q\ in our 
basis. On the lattice, this operator mixes under renormalization with the other four 
operators as follows 

r 5 
Qi = Z11 Qi + ZuQ 

i=2 

where Z\\ is a multiplicative logarithmically divergent renormalization constant; it 
depends on the coupling and a//. The mixing coefficients Z\i (with i = 2, ... ,5) are 
finite; they only depend on the lattice coupling <?o( a ) 

The renormalization of the parity-conserving parts of the operators O^ 2 and O^ 2 
is related to that of the operators Q2 and Q3; the correspondence is given by 



3 7 /2 -» Q 2 , 
0l' 2 -> -2Q 3 - 

The renormalized operators are defined as: 

Q2 = Z22Q2 + Z23Q3, 

Q3 = Z32Q2 + (8) 

where 2L- (with i,j = 2,3) are logarithmically divergent renormalization constants 
which depend on the coupling and a\x. The above mixing matrix is not peculiar to 
the lattice regularization, but also occurs in the continuum. The breaking of chiral 
symmetry by the Wilson action requires the additional subtractions: 

Qi=Qi+ E Z »Qh » = 2,3. (9) 

.7=1,4,5 
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where the ZijS are finite coefficients which only depend on <?o(a). 

Finally, the operator O^J 2 corresponds to the operator Q\ of eqs. (||). Thus, its 
renormalization constants, ^-parameter, etc. are identical to those of 0^ s=2 . The 
results for all the renormalization constants Zij (computed with the NPM at several 
renormalization scales fj, at (3 = 6.0) can be found in [O]. 



4 B- parameters 

In the AS = 2 case, the S-parameter is defined as 



K{H (K°\6*s=*\K°) VSA ' 1 > 

where the matrix element in the VSA is given by 

(K°\d* s = 2 \K°) VSA = 2 (l + £) Z 2 A \(K°\A,\0)\ 2 = 2 (l + £) r^/i , (11) 

with the (finite) renormalization constant of the lattice axial current = ^7 . 
Since (K°\O as=2 \K°)vsa is given in terms of physical quantities (txik and fx), B(fi) 
runs with the scale \i exactly like the corresponding renormalized operator 0(/x). An- 
other definition of the VSA matrix element has recently been proposed in ref. ||. It 
consists in vacuum-saturating each operator which enters in the subtraction of eq. (0). 
In this way, statistical fluctuations are reduced in the ratio of correlation functions 
used to extract Bk (defined as R AS=2 in eqs. ( |l5|) below). As pointed out in ref. |27|| , 
however, the definition used in ref. Q spoils the good scaling properties of Bk{/J>)- 



Thus, we insist on retaining the standard definition of eq. (11), at the price of having 
larger statistical errors. 

For the AI = 3/2 transitions, the i?-parameters are defined by: 

(n+\0] ,2 \K+) VSA 
(Tr+lOg 7 \K+)vsa 

As discussed in sec. ||, B^ 2 ([i) = Bk{h) in the limit of degenerate quark masses. The 
VSAs for the above matrix elements depend on two different contributions of the form: 

(7r + \d 3 7 /2 \K + ) VSA = ^Z 2 P (7r + \P\0)(0\P\K+)-Z 2 A (TT+\A ll \0)(0\A tl \K + ), 

(tt+\6 3 8 /2 \K+) V sa = 2Z 2 (tt + \P\0)(0\P\K + )-^(tt + \A i1 \0}(0\A i ,\K + }, (13) 
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where Zp is the renormalization constant of the lattice pseudoscalar density P = ^75^ 
(renormalized at the same scale \x in the RI scheme). Since we work with degenerate 
quark masses, we have left the flavour content of the operators and P unspeci- 
fied; they are meant to have whatever flavour is required by the hadronic states of 
their matrix elements (i.e. P n , Pk and similarly for An)- Contrary to the AS = 2 
case, the leading terms of the above VSA matrix elements are /^-dependent quantities, 
with an anomalous dimension equal to twice the anomalous dimension of the pseu- 
doscalar density P. Thus the B-parameters do not scale in fj, like the matrix elements 
of the corresponding operators. The last terms on the r.h.s. of eqs. (^) vanish in 
the chiral limit. Consequently, since we are ultimately interested in passing from the 
{ir + \6 3/2 \K + ) matrix elements to the <7r7r|0 3/2 | K) ones using soft pion theorems, fol- 
lowing ref. [16], we have dropped the last terms on the r.h.s. of eqs. In order 



to extract the ^-parameters, we need to compute the following two- and three-point 
correlation functions: 

G P (t x ,p) = Y / (P(x)P(0))e-^ , G A (t x ,p)=Y / (Mx)P(0))e-^, (14) 

X X 

G d (t x ,t y ;p,q) = Yl(P(y)6(0)P(x))e-^e^, 

x,y 

where x = (x,t x ),y = (y,t y ) and O stands for any four-fermion operator of interest. 
As stated above, all correlation functions have been evaluated with degenerate quark 
masses and therefore only the "eight-diagrams" contribute to Gq. By forming suitable 
ratios of the above correlations, and looking at their asymptotic behaviour at large 
time separations, we isolate the desired operator matrix elements. In particular the 
ratios: 



R 



AS=2 



1 G 6AS =2 (K°(p)\6 AS = 2 \K (p)) 



Z\ GpGp Zl\(0\P\K°)\< 



3/2 N c G or N c (vr+|d 7 3/2 |A+) 

7 2Z 2 P GpG P ^ 2 Z 2 P (n+\P\Q)(0\P\K+y 1 ' 



R 



.3/2 1 G 6 8 3/2 1 (n + \6l /2 \K- 



2ZpG P Gp 2Z^+|P|0)(0|P|A+) 



give, up to computable factors, the S-parameters of interest. For comparison, we have 
obtained results with the operator renormalized not only with the NPM, but also with 
SPT and BPT. 



5 Numerical results 

Our simulation has been performed at (3 = 6.0 with the tree-level Clover action in 
the quenched approximation. Quark masses have been fixed at k = 0.1440,0.1432 and 
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0.1425. The renormalization constants have been obtained from quark correlation func- 
tions, in the Landau gauge, on a 16 3 x 32 lattice, with 100 configurations. The hadronic 
matrix elements have been computed on an 18 3 x 64 lattice with 460 configurations. 
Details on the choice of time intervals, spatial momenta and related technicalities are to 
be found in ref. IfT^] . Statistical errors have been estimated with the jacknife method, 
by decimating 46 configurations at a time. We have neglected the statistical errors of 
the renormalization constants, quoting only those of the hadronic matrix elements. In 
the above ratios, we also need the (finite) axial-current renormalization constant Za 
and the o^-dependent renormalization constant Zp of the pseudoscalar density. De- 
pending on the method of renormalization of the four-fermion operator (NPM, SPT or 
BPT), we have used the corresponding estimate of Za and Zp, obtained with the same 
method of calculation. Although Za should not depend ^ on afi, slight variations of its 
NPM estimate, arising from systematic effects, partially cancel analogous variations of 

R AS=2 

giving much more stable results. The NPM estimates for Zp and Za used in 



the present work are those of ref. [29| 



In order to extract the £>-parameters from the ratios of eqs. fllq), we follow the 
procedure of ref. 12], fitting them with the function 



where 



R = a + (3X + 7 y , (16) 



x = gg# « J!±„,. r = ^x (17) 



3G P Gp 3Zl|(0|P|K»)| 2 



m 



K 



(the large time asymptotic behaviour of X is also shown above). The results of the fits 
are collected in tab. [I], at eight different values of the renormalization scale /j, 2 a 2 . We 
also show results obtained with the renormalization constants evaluated in SPT and in 
BPT. In the latter case, we use the recipe a b ° ost = a s /(^Tr(Un)), where (Ua) stands 
for the average plaquette. At (3 = 6.0, we have used a h ° ost = 1.68/4-7T. 

We first examine the results for R^ s=2 , from which Bk can be extracted. As 
in [p4|]-[13], it can be seen that a and (3 (which are lattice artifacts) are compatible 



with zero within at most 2<r, unlike their SPT and BPT values. Thus (within our 
statistical accuracy and up to terms of 0(g 2 a)), the correct chiral behaviour of the 
matrix elements is restored using the NPM. We point out the stability of 7 as a function 
of fi. Assuming a and (3 to be zero, the -B^-parameter is then given by: 

B K = 7- (18) 

In our original analysis JTT| , p^ , we have erroneously considered the mixing of 
the AS = 2 operator with only the three other operators which appear in one-loop 
perturbation theory M. It has been stressed in g [H| [14|, Q that, non-perturbatively, 



2 For the tree-level Clover action, the best available non-perturbative estimate for Za is obtained from 
lattice axial WPs; at (3 = 6.0 this is Z A = 1.11(2) |1. 
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Operator 


H 2 a 2 


a 


P 


7 




0.31 


0.027 ±0.014 


0.17 ±0.16 


0.68 ± 0.12 




0.62 


-0.018 ±0.014 


0.28 ± 0.17 


0.67 ±0.12 




0.96 


-0.014 ±0.014 


0.27 ±0.16 


0.66 ±0.11 




1.27 


-0.010 ±0.013 


0.26 ± 0.16 


0.66 ±0.11 




1.39 


-0.004 ±0.013 


0.23 ± 0.16 


0.66 ± 0.11 


R AS=2 


1.85 


-0.005 ±0.013 


0.25 ± 0.16 


0.66 ± 0.12 




2.46 


0.002 ± 0.013 


0.25 ± 0.16 


0.67 ± 0.12 




4.01 


U.Ulz ± U.U12 


0.25 ± 0.1b 


n co i nil 

O.o8 ± 0.12 




SPT 


—0.069 ± 0.013 


0.17 ± 0.16 


0.65 ± 0.12 




BPT 


—0.058 ± 0.013 


0.18 ± 0.1b 


O.DD t 0.12 




0.31 


1.70 ±0.16 


-6.4 ± 1.3 


6.3 ±0.7 




0.62 


0.81 ± 0.07 


—0.8 ± 0.6 


3.5 ± 0.4 




0.96 


0.69 ± 0.04 


—0.1 ± 0.4 


2.6 ± 0.3 


R 3 7 ' 2 


1.27 


0.70 ± 0.03 


0.1 ± 0.4 


2.3 ± 0.3 


1.39 


0.70 ±0.03 


0.04 ± 0.34 


2.2 ± 0.3 




1.85 


0.69 ±0.03 


0.15 ± 0.28 


1.94 ±0.23 




2.47 


0.69 ± 0.02 


0.22 ± 0.25 


1.77 ± 0.21 




4.01 


0.73 ±0.02 


0.31 ± 0.22 


1.62 ±0.20 




SPT 


0.68 ±0.02 


0.45 ± 0.17 


1.14 ±0.14 




BPT 


0.48 ±0.02 


0.41 ± 0.22 


1.58 ± 0.17 




0.31 


1.20 ±0.06 


-1.0 ±0.4 


0.66 ± 0.20 




0.62 


1.08 ±0.04 


-0.5 ±0.3 


0.56 ± 0.15 




0.96 


1.04 ±0.03 


-0.3 ±0.3 


0.53 ± 0.14 


R S/2 


1.27 


1.00 ±0.03 


-0.06 ±0.25 


0.52 ± 0.13 


1.39 


1.00 ±0.03 


-0.05 ±0.25 


0.52 ± 0.13 




1.85 


0.99 ±0.03 


0.04 ±0.24 


0.51 ± 0.13 




2.46 


0.99 ±0.02 


0.09 ± 0.24 


0.51 ± 0.13 




4.01 


0.98 ±0.02 


0.15 ± 0.23 


0.51 ± 0.13 




SPT 


0.81 ±0.02 


0.36 ± 0.17 


0.42 ± 0.09 




BPT 


0.75 ±0.02 


0.29 ± 0.16 


0.39 ± 0.09 



Table 1: Values of the fit parameters for R AS 2 , R 7 and R s , from the NPM (at several 
renormalization scales), SPT and BPT. 



qAS=2 a i go m j xes w ith a fourth operator. The results of ref. [13] and the present work 
take into account the complete non-perturbative mixing. 

We now turn to Rj , R^ 2 , from which B^ 2 and B^ 2 can be extracted. Note the 
stability of a as a function of a^i and the compatibility of (3 with zero. In this case, for 
both operators, the S-parameters are simply given by 

5 3/2 = a. (19) 

Comparing the NPM results to those obtained in SPT and BPT, we find agreement 
between the NPM and SPT values of B^ 2 , whereas its BPT value is incompatible with 

3 11 

the other two. The B s result depends on the method used for its renormalization. 

6 Physics results and conclusions 

The S-parameters which can be extracted from the results given in table [l] have been 
obtained in the RI scheme, since the operators have been renormalized with the NPM. 
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As previously stated, it is customary to express all results in the MS scheme (with 
NDR dimensional regularization) . The finite matching between the RI and MS-NDR 
renormalization schemes can be done in continuum perturbation theory, by computing 
the operator matrix element in the same gauge and on the same external quark states as 
those used for the non-perturbative calculation of the lattice renormalization constants. 



At NLO this matching has been obtained in [22|. We have used it in order to convert 



our NPM results from the RI to the MS-NDR scheme. The SPT and BPT cases are less 



straightforward: we have started from the renormalization constants of ref. |3C], which 
relate the lattice operators to those renormalized in the MS-DRED scheme. To these 
constants we have subtracted the contributions due to the use of the DRED scheme 
and added those corresponding to the RI one. In this way we have obtained, for both 
the SPT and BPT cases, the renormalization constants in the RI scheme, from which 
the results of tab. [l] are computed. Alternatively, we have also calculated and added 
the necessary MS-NDR contributions in order to obtain the SPT and BPT results in 



this scheme. The latter calculation has also been performed in 14] for the Wilson 
action. At NLO j22|, the result for B K is 

BWM = ( 1 - ^Arf) Bf(»), (20) 



where 

Ar^ s = 14/3-81n2. (21) 

The matching relation for the renormalized operators O, (i = 7, 8) at NLO is also 



given in |2_£ 



where 



Ar Ms_ ( § + §1*2 -2-2M2 \ 

Ar ^' " 2-21n2 + ' {2S) 



Hence, for the B-parameters B^ 2 and B^ 2 we have: 
( B 3 7 /2 \ _ 1 / 1 - ^?Arff _AT c ^Ar|S 

\ R 3 / 2 / / nfuK T7c\2 1 1 &sW A r MS -I aAM} a„MS I I R 3/2 

V B 8 /MS (l + ^-Arf S ) \ ~A^T^ Ar 87 1 --5T Ar 88 / V #8 / rj 

(24) 




with 



Arf s = 16/3. (25) 



The uncertainty due to the choice of a s ((i) is about ±0.03. We find a further uncer- 
tainty of ±0.03 when varying the number of active flavours from to 4. 

In order to compare our result for Bk with other theoretical predictions of the same 
quantity, it is useful to convert it to the RGI quantity, Bk, by multiplying it by the 
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Wilson coefficient. At NLO, Bk is given by 

a.(jj,) / 7 (1) /3o-7 (0) A' 



B K = a s (n) 



-7 (0) /2/3o 

4tt I 2(31 



B™ s (ri, (26) 



where /?o,i and 7( 0,1 ) are the leading and next-to-leading coefficients of the /3-function 
and anomalous dimension. /3q 1 and 7^) are universal whereas 7W depends on the reg- 
ularization and renormalization schemes. The explicit expressions of these quantities 
can be found for example in |^]. This estimate of Bk is also regularization-scheme 
independent, up to next-to-NLO order terms. 

In tab. |2] we collect all our results for the _B-parameters, obtained by computing 
their renormalization with the NPM and in BPT, at the reference scale [i ~ 2 GeV in 
the MS scheme. For Bk we also give the RGI value Bk- Moreover, we compare our 
results to the most recent ones obtained with the Wilson action at the same coupling 
and with operator renormalizations carried out perturbatively in BPT J[4|]. 

Two points concerning this comparison deserve some attention here. The first is 
that in ref. JIJ|] the results were obtained in 1-loop BPT by renormalizing the lattice 
operators directly in the MS-NDR scheme, rather than doing the renormalization in 



the RI scheme and then matching to MS-NDR with eq. (24), as we have to do with 
the NPM. In perturbation theory, the two procedures differ by O(a^) terms, which 

introduce a systematic difference of about 12%-15% in the results for B 7 and B 8 . 
For this reason, in order to make the comparison meaningful, our BPT estimates in 
the MS-NDR scheme, given in Table ^, were obtained using the same procedure as in 
|14| . The second point is that in ref. []14| , a more complicated version of BPT than in 
this work has been implemented, which involves tadpole resummations and a choice 
of "optimal" scale q* for the boosted coupling. Since the recipe for q* is not unique, 



results with several choices of q* have been listed in [14], and can be compared with 
our results. Finally, we point out that a direct comparison of our Bk value to those of 
j|] is not possible, because the latter have been obtained at different lattice couplings. 

Our results from NPM and BPT for Bk and Bk are in perfect agreement. With a 

larger statistical error, our Bk value also agrees with those of refs. [H, |31j| . We find, 

3/2 3/2 

instead, a discrepancy between our NPM and BPT estimates of B? and Bg . Our 
values obtained with BPT are fully compatible to those of JTJ] (at least for one q* value), 
where the Wilson action was used. The NPM estimate, instead, is in disagreement 
with any value obtained in BPT (either with the Wilson or the Clover action and 
for several boosting variants). This indicates that the difference between our NPM 
estimate and that of ref. |l4| is due to the NPM used in the former result, rather than 
the implementation of different actions (Clover and Wilson respectively). The increase 

3 /2 

in the value of B 8 , obtained with the NPM, is of great phenomenological interest, 
since it may induce a significant decrease of the ratio e'/e. 

In conclusion, the non-perturbative renormalization of the four-fermion operators, 
either with the NPM or with WI's, strongly improves the reliability of lattice com- 
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Quantity 


Method 


Result 


Ref. 


Bk 


NPM 
BPT 
BPT q* = 1/a 


0.66(11) 
0.65(11) 
0.74(4) 


this work 
this work 

R 


Bk 


NPM 
BPT 


0.93(16) 
0.92(16) 


this work 
this work 




NPM 

BPT 
BPT q* = 1/a 
BPT q* =n/a 


0.72(5) 
0.58(2) 
0.58(2) 
0.65(2) 


this 
this 


w 
w 

u 

11 


ork 
ork 


R 3/2 


NPM 

BPT 
BPT q* = 1/a 
BPT q* = vr/a 


1.03(3) 
0.83(2) 
0.81(3) 
0.84(3) 


this 
this 


work 
work 

la 



Table 2: B-parameters for AS = 2 and Al = 3/2 operators at the renormalization scale 
fi = a^ 1 ~ 2 GeV. All results are in the MS scheme. 



putations of l?-parameters with Wilson- like fermions. Perturbative mixing of lattice 
operators, carried out at lowest order, fails to reproduce the expected chiral behaviour 
of the matrix elements (at present day couplings). A good chiral behaviour of 0^3=2 

3 /2 

is restored when the operators are renormalized non-perturbatively. Moreover, B 7 

3 /2 

and B 8 obtained with the NPM disagree by as much as 20% from the ones obtained 
in BPT. 
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